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Abstract 

We discuss the ability of the generator coordinate method (GCM) to select 
collective states in microscopic calculations. The model studied is a single- 
j shell with hamiltonian containing the quadrupole-quadrupole interaction. 
Quadrupole collective excitations are constructed by means of the quadrupole 
single-particle operator. Lowest collective bands for j=31/2 and particle num- 
bers A^=4,6,8,10,12, and 14 are found. For lower values of j, exact solutions 
are obtained and compared with the GCM results. 

PACS Numbers 21.10Re, 21.60Ev 



1 



Typeset using REVTgX 



I. INTRODUCTION 



A microscopic description of nuclear collective modes is still a very important and in- 
teresting problem studied in the theory of nuclear structure. In a frame of the shell model 
such a description is difficult to obtain because of a very large size of the basis required to 
include all necessary configurations. For example, only very recently excellent shell-model 
results for rotational states in ^^Cr have been obtained Although many phenomeno- 
logical approaches can, and have been constructed, their microscopic derivation from the 
shell model is often not entirely accomplished. It is therefore important to analyze and 
study approaches which are able to provide a simultaneous description of collective and 
single-particle phenomena. 

The generator coordinate method (GCM) constitutes an effective method to deal 
with collective degrees of freedom while starting from a pure many-fermion description. 
Moreover, it may also yield an exact quantum mechanical formalism, and only the choice of 
generating functions and a generator coordinate decides whether we can find exact solutions, 
collective as well as noncollective, to the actual problem. The configuration mixing feature 
which is by construction built in into the GCM provides us with the ideal tool to study 
microscopic foundations of collective models. 

In the present study, we apply the GCM to a microscopic model of a single-j shell filled 
with even number of identical nucleons interacting via pairing-plus-quadrupole hamiltonian 
1^. Although the model is very simple, it contains two main ingredients of the "real" nuclear 
structure theories, i.e., the quadrupole deformation and the pairing correlations, and at the 
same time for low values of j it is relatively easy to solve. On the other hand, the GCM is 
well suited to study this model for large values of j where the exact solutions are inaccessible. 

In the past the single-j shell model has been studied by many authors. Mulhall and Sips 
found that even for 4 particles occupying single-j shell with strong quadrupole-quadrupole 
force the collective effects are important and rotational structures are present in the spectra. 
Deformability, i.e., a competition between the quadrupole force and the pairing, was studied 



in a series of papers by Baranger and Kumar 13,0]. Arima [10| discussed excited rotational 



level structures {K=2 bands) appearing in the single-j shell. Friedman and Kelson fTL 



analyzed collective spectra as depending on particle number and j. The present paper aims 
at using the GCM as a filter which would select collective structures among the complete 
shell-model spectrum. Such a study may give us more confidence in using the GCM in more 
realistic cases where the quality of approximations involved is difficult to analyze. 

The paper is organized in the following way: In Sec. |T|we fix notation by briefiy reviewing 
the pairing-plus-quadrupole model in a single-j shell and then in Sec. [II A| we present the 
available exact solutions. The key point of the GCM is the construction of the generating 
states, which we present in Sec. B| by invoking the single particle coherent excitation model 
(SCEM) of Dobaczewski and Rohozihski In this model the quadrupole excitations 

are build using the single-particle quadrupole operator. In the present study we restore 
the broken particle-number and angular-momentum symmetries exactly. This is done in 
a consistent GCM framework by using the gauge angle and the Euler angles as generator 
coordinates. The particle-number projection [|14| , |T5| is based on the Fomenko method |ll6|JT7[| 
in which integrals are replaced by finite sums. Sec. [Ill A| , and a similar method is also used 
to perform the angular momentum projection [[T8[,|l9|, Sec. [Ill B[. A description of the GCM 
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calculations in the intrinsic frame of reference is presented in Sec. In Sec. [IV A| we 
analyze the GCM results for j=15/2 and A^=8 particles and compare them with the exact 
calculations. In Sec. |1V B| we present the GCM spectra for j=31/2 and N=4, and for particle 
numbers between 6 and 14, where the exact solutions are not available. 



II. SINGLE-J SHELL 

We consider the (2j+l)-fold degenerate single shell of angular momentum j filled with 
even number of identical particles, which without the interaction is assumed to be at zero 
energy. The hamiltonian is composed of the pairing-plus-quadrupole interaction, 

H = -GP+P-xQ-Q, (2.1) 
where is the pair transfer operator and Q is the quadrupole moment operator, 

P+ =Y,{jmjm'ma+a+, , (2.2a) 

mm' 

Qt = IZ(i"^^'"^'|2/i)a+a^/ , (2.2b) 

mm' 

while G and x ^i^re pairing and quadrupole coupling constants, respectively. Hamiltonian 
( p.l| ) describes basic collective correlations between nucleons and it has been used by 



many authors §-11,20,21 



In the present paper we aim at solving the pairing-plus-quadrupole model in a large 
single-j shell. Dimensions of the many-fermion space increase rapidly with j and the exact 
solutions become inaccessible very fast. For example, in the M and J representations the 
maximum dimensions for j=15/2 are 526 and 35, respectively, while for j=31/2 they are as 
large as 8 908 546 and 200 691. The generator coordinate method, which we apply in Sec. 

to this problem, is able to describe low lying collective states without using very large 
matrices. Here we present solutions for the values of particle number N and single-particle 
angular momentum j for which in the M representation the dimensions of the many-particle 
space are not larger than 5000. 



A. Exact diagonalization 

In Fig. [| we show the yrast bands for A^=8 particles in the j=15/2 shell obtained with 
several values of the coupling constants G and Xi see Eq. |2.1| . Using the parametrizations 
G={l—x) X 1 MeV and x=xx 1 MeV we show a transition from the pairing (seniority) limit at 
a;=0 to a large-deformation limit at x=l. It can be seen that for x<l the spectra are strongly 
influenced by the flxed-seniority structures which at x=0 appear as degenerate multiplets. 
Since in this paper we want to discuss situations corresponding to collective quadrupole 
excitations, in the following we present only the results for the pure quadrupole-quadrupole 
force, i.e., we set G=0 and x=l MeV. 

In Fig. |](a) we present the yrast bands for A^=4 particles interacting with the quadrupole- 
quadrupole interaction in shells corresponding to j between 15/2 and 31/2. At low spins, 
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regular rotational-like bands are obtained with moments of inertia jr=3^^/i?2+ varying be- 
tween 51 ^^/MeV for j=15/2 and 333 ^^/MeV for j=31/2, and the corresponding E4+/E2+ 
ratios varying between 3.23 and 3.32 {Ej+ denotes the excitation energy above the ground 
state). 

At higher spins the ground-state bands are crossed by bands built on high-i^" particle- 
hole excitations, which can easily be identified with oblate structures. In the Nilsson 
single-particle diagram presented in Fig. |^ for j=15/2, we see that at the oblate side the 
four particles occupy the ^2=15/2 and 13/2 orbitals. Therefore the lowest single-particle 
excitation corresponds to promoting a particle from Q=13/2 to f2=ll/2 which gives the 
7^^=13/2+11/2=12 excitation and leads to a band head at 1=12. Similarly, for larger values 
of j one obtains band heads at larger even values of spin. 

For A^=6 particles. Fig. |^(b), we see the bands which correspond to 2p-2h oblate configu- 
rations. For example, for j=21/2 the p-h band corresponds to exciting the ^2=17/2 particle 
to the f2=15/2 orbital {K=16) and then the second p-h excitation aligns the f2=19/2 parti- 
cle to the f2=13/2 orbital leading to i^"=32 and giving the band head at /=32. We see that 
for small shell-filling factors A^/(2j+l) the high-ii' bands appear at very low energies, while 
they cease to be yrast for larger shell-filling factors. For example, for N=8 and 10, Figs. H(c) 
and (d), no yrast high-i^ bands are seen and only yrast-yrare band interaction is visible at 
7=12 for the lowest j=19/2 and A^=8 bands. In general, the yrast bands for half-filled shells 
present regular collective structures, although they do not have pure rotational character, 
for example, i?4+/i?2+=l-92 for the j=A^— 1/2=15/2 band. This is due to the 7-instability 
of the corresponding potential energy surface (PES), see Fig. ^(a). 

The complete spectrum of 526 states for A^=8 particles in the j=15/2 shell is presented 
in Fig. ^. The lines connect states with the largest E2 reduced matrix elements of the 
quadrupole operator. Only the stretched AJ=2 matrix elements are considered. The lowest 
O'^ state is connected to the 2+ state with which it has the largest E2 matrix element. Then 
the first excited 0"*" state is connected to one of the remaining 2"*" states, and the procedure 
is continued for still higher O"*" states. After that all remaining 2"*" states are considered to be 
the band heads. In the next step the 2+ states are considered by order of increasing energy 
and connected to the 4"*" states. In this way all even spins are connected into bands and 
then the same procedure is repeated for odd spin values. 

Several regular collective bands can be seen in Fig. ||. These will be used to discuss the 
ability of the GCM to select collective structures in the sea of all possible eigenstates. Sec. 
[IV A| . Here we only discuss a possible interpretation of lowest states in terms of the collective 
model. In Table | we present the reduced E2 matrix elements between five lowest O"*" and 
2"*" states. It can be seen that these matrix elements obey a strong selection rule related 
to the special symmetry of the half-filled shell. Indeed, the quadrupole operator ( p.2b| ) is 
odd with respect to changing the creation operator into the annihilation operator and vice 
versa. Therefore, the quadrupole-quadrupole hamiltonian is invariant with respect to such a 
change and the spectra of systems containing N and 2j+l— particles are strictly identical. 
However, for a half-filled shell A^ is equal to 2j+l— A^ and the above symmetry allows to 
attribute a new dychotomic quantum number to every eigenstate, while the quadrupole 
operator may only connect states with these quantum numbers being different. 

An approximate interpretation of the lowest states can be done in terms of the collective 
Jean-Wilets model p2[, which describes collective states for the 7-unstable PES. In this 
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model the lowest-order quadrupole transition operator changes the seniority quantum num- 
ber by one and therefore the parity of the seniority can be associated with the dychotomic 
quantum number discussed above. In this way the O^*" and states can be identified with 
seniority zero and three and the three lowest 2"*" states with seniority one, two, and four. In 
this scheme the states 0^ and can be interpreted as mixtures of the seniority six states 
with the seniority zero (3 vibrations, and the states 2j and 2^ as mixtures of seniority five 
and seven with seniority one f3 vibrations. This interpretation is however not manifested in 
the energy spectra, where the characteristic Jean-Wilets multiplets are strongly split. One 
can also note that if the identification of the (3 vibrations is correct, the corresponding energy 
is rather high as compared to the energy of the motion in the 7 direction. This high energy 
may also suggest that the collective 7 vibrations are here directly mixed to non- collective 
states and that the f3 vibrations are absent in the model. 

The first excited band can be interpreted as the 7 vibration strongly coupled to rotations 
due to the 7-instability of the corresponding PES, Fig. ^(a). At spin /=8 it is perturbed 
by a coupling to the lowest noncollective excitation promoting the Q=7/2 particle to the 
Q=9/2 state (if one considers the oblate side of the Nilsson diagram. Fig. ^). However, one 
is here unable to distinguish between the oblate and prolate structures, because the prolate 
side suggests an excitation of the ^2=9/2 particle to the Q=7/2 state, which gives the same 
value of spin 7=8. 

For j=31/2 the exact solutions are available for N=4 only. Fig. The collective bands 
are identified by the same procedure as described above using the reduced matrix elements 
of the quadrupole operator. At low spins one can see several regular bands which can be 
fairly well interpreted in terms of the 7 vibrations coupled to an axial rotor with K being 
approximately a good quantum number. In this interpretation, the ground-state band is 
the K=0 and n^=0 band and the first excited band has K=2 and n^=l. Then come two 
n-y=2 bands with K=4 and K=0 followed by two n^=3 bands with K=6 and K=2. One 
may even identify three n^=4 bands with K=8, 4, and 0. Again, in this interpretation there 
is no room for the (3 vibrations at low energies ||23|. Table |I| shows the E2 reduced matrix 



elements between the lowest 0"'" and 2+ states and confirms the above interpretation. A clear 
decrease of the matrix elements with the difference in the number of 7-vibrational quanta is 
manifest, however, the strict |An^|<l selection rule is not present. 

Similarly as for the yrast bands. Fig. 0, one may easily recognize in the spectrum of Fig. 
^ the bands built on the oblate particle-hole excitations. Together with the excitation of the 
f2=29/2 particle to the Q=27/2 orbital, which becomes yrast at 7=28, one can also identify 
the excitations of the same particle to the Q=25/2, 23/2, 21/2, and 19/2 orbitals. One may 
clearly see the effect of decreasing E2 probability when the band terminates; the band heads 
are then coupled stronger to adjacent bands then to the next members of the given band. 



B. The SCEM model in a single-j shell 

The single-particle coherent excitation model (SCEM) proposed by Dobaczewski and 
Rohozihski ||T2| , p!3| is based on the assumption that the quadrupole excitations can be con- 
structed by acting with the quadrupole fermion operator 7^+ on a reference state of the 
form 
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|ref) = (^+)^/' |0) , (2.3) 

i.e., 

|^scEM) = (i^+)'^" |ref), (2.4) 

where creates a monopole collective pair of particles and is a single-particle 
quadrupole excitation operator. Since only one operator of a given mutipolarity can be 
constructed in the single-j shell, here S~^=P^ and F^=Q. In the SCEM, for Np=0 we 
obtain one state of 1=0 and for Nf=1 one state of 1=2. If for every value of Np we orthog- 
onalize the SCEM states with respect to those obtained for Np—l, we obtain the numbers 



of states listed in Table |T| for j=15/2 and N=8. 



It can be seen that for low values of Np the numbers of states of a given / exactly 



correspond to those of the five- dimensional harmonic oscillator |24| with Np phonons, or to 
the numbers of d'^s excitations required to create the interacting-boson-model states ||2^ 
from the s-boson condensate {s'^)^^ . However, SCEM being a true fermion model, there 
are departures from this rule appearing as soon as Np reaches half of the shell's maximum 



particle number j+1/2 (i.e., in Table |nT| from Np=A on), which reflect the effect of the Pauli 
correlations. It is obvious that states of the angular momentum I can only be obtained for 
Np>I/2 and therefore the number of excitations Np plays the role of the cut-off in spin. 

The generating function of the GCM is here constructed as a condensate of the 
excitations acting on a condensate of the pairs 

I tot2) = exp(t2 ■ F+) exp(to^+) | ) . (2.5) 

The parameters to and t2 are a scalar and a quadrupole numerical tensors, respectively, and 
are used as the generator coordinates. We use the notation of the scalar product defined as 
t2-F^=J2tit*fj,^ll and assume that t*^={—l)'^t2-fj, and to=^o- 

For the GCM state ( p.5|) neither the particle number nor the angular momentum is 
a good quantum number. It can, however, be decomposed into the SCEM states ( p.4|) 
which, with a proper angular momentum coupling of the excitation operators F~^, have 
good angular momentum and particle number. The average value of the particle number in 
the generating function ( |2.5| ) increases with the to parameter. Similarly, the average number 
of the excitations F^, and hence the average value of the angular momentum, increases with 
the magnitude of the quadrupole tensor 



/5' = El^2,r (2.6) 

It is convenient to introduce three matrices 

S^m' = {jmjm'm = {-iy-"^^^6m,-m' , (2.7a) 

V2j + 1 

F:l> = ijmj-m'Mi-iy^"^', (2.7b) 

^mm' ~ E^2At-^mm' ' (2-7c) 

where S is antisymmetric and F is hermitian, in terms of which the basic SCEM building 
blocks read 
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— ^ Smm'O'm^m' i (2.8a) 
mm' 

= E Km'<<^n^' , (2.8b) 
mm' 

h ■ F'^ = E Fmm'O^m'^m' ■ (2.8c) 
mm' 

Now, the generating function ( ^^.51 ) can be represented as 

|tot2) = exp(toe*-^''^+e-*-^^) |0) , (2.9) 



and further, using the nonunitary Bogohubov transformation ||2^, as a Thouless [g^ state: 



where 



\toh) = exp(i C^m'«') |0) , (2.10) 



C^^m' = 2toT.ie^)^nSL'{e^ )n'm' ■ (2-11) 



In this way all GCM kernels can be easily calculated using the standard expressions resulting 
from the Wick theorem |2^J^. On the other hand, the product of matrices in ( p.ll|) can be 
rapidly obtained by using the F matrix in its eigen-reference-frame, which corresponds to 
transforming the tensor t2 to the intrinsic system defined by conditions 

^21 = , ^22 = ^22 ' (2.12) 

and then diagonalizing the mean field F in the intrinsic frame. In cases when the 5*"^ pair 
has the seniority form, i.e., contains pairs of time- reversed states with equal amplitudes as 
is in the single-j shell ( p.7aj ), we may further simplify expression ( p.ll| ) by using the fact 



that F is a time-even matrix, i.e., 



^mm' — \. 1 ^^"^^^m-m' ■ (2-13) 



III. RESTORATION OF BROKEN SYMMETRIES 

It is well known that the GCM for the gauge angle and the Euler angles used as generator 
coordinates is equivalent to the exact particle-number and angular-momentum projections 
P(]|,|^J5|]. In the present study we use this property to restore those broken symmetries 
exactly, and then the GCM equation is solved in the intrinsic frame of reference. 
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A. Particle-number projection 

A standard way [|14|,|15|,^^ to perform the particle-number projection is to introduce a 
gauge angle and to integrate over its (0, 2it) domain. Using the generating function in the 
form (|2.10|) , the projected state is 

\Nt,) ^ pM\tot2) = ^ fd<P e-^*exp(i ^ C^^>e'^'^aUi)\0) . (3.1) 

^ mm' 

Such an integration leaves only the components of the wave function with the particle number 
N. This is based on the following orthogonality relation valid for any integer N and N' 

^ rdcl> e^(^'-^)'^ = 5^^, . (3.2) 



27rJo 

When calculating matrix elements of any arbitrary particle-number-conserving operator O 
between the Thouless states, the integration over the gauge angle can be done only for one 
of those states 

{Nt2\d\NQ = {tohlOpNKt',) = ^ rd<f) e-''''^{tot2\d\e'%,t',) , (3.3) 

ZTT Jo 

which amounts to multiplying the Thouless matrix C (or, in other words, the coefficient t'o) 
of the ket state by the factor e^**^ and integrating over 0. 

When there is an upper limit A^max on the particle number of the system, as is the case 
in the single-j shell, the integral ( p. 3] ) can be discretized as proposed by Fomenko |]TB|,|T7 . 
Instead of Eq. ( p.2|) one can write 

r 1 for iV' = iV 

1^1 ^0 ^ ' l-exp[.(iV- iV)(iV +1)A0] ^ 

(3.4) 

Taking A0=27r/(A^max+l), the numerator of the expression for N'^N vanishes for all 
and N', whereas the denominator does not vanish for any and A^' not larger than A^max, 
i.e., 

1 A^max 

' . E e^(^'-^)"^<^ = 5^;v' for AT, AT' < AT^^x • (3.5) 

Amax + -L n=0 

For systems with even particle numbers we may repeat a similar discretization using only 
|A^max points and A'0=vr/(|Amax+l), which gives 

1 

1^ \, E e*(^'-^)"^'^ = 5^^, for even A, A' < A^^x • (3.6) 

2 max + J- n=0 

Furthermore, assuming that |Amax is also even, or proceeding as if Amax was larger by 2 
when ^Amax is odd, the left hand side of Eq. ( p.6| ) can be rewritten in a form of the sum 
over n from -jN^^^ to ^A^ax, and finally 



8 



-N +12 

2 * max J- L i 



^ -^max 



^j(Af'-Ar)nA'0 



n=l 



Jtvtv/ for even A^, A^' < A^^ 



(3.7) 



Since the unprojected GCM states are time-even, the matrix elements of time-even operators 
are real and then the prescription for the projected matrix elements takes the following form 



{Nt2\d\Nt'2) 



1 



{tot2\d\t'QQ + ^ 



^ ^ may 

E 

n=l 



-iNnA'cf) 



(3.8) 



where only jNj^ax integration points are required. 

The standard formulas resulting from the Wick p9|j5[l and Thouless [^,0 theorems ap- 
plied to evaluate (to^2|0|e^*"''^ '^t'^, t'2) matrix elements contain the products of matrices C^C , 
cf. Eq. ( p.lO| ). Since using the projection method requires many such matrix multiplications 
a very efficient algorithm can be obtained by transforming C^C to the basis in which it is 
diagonal (Appendix 0). Then the matrix multiplications reduce to simple sums and prod- 
ucts of the corresponding eigenvalues. At the same time the sign problem of the GCM 



kernels is avoided by considering only one half of pairwise degenerate eigenvalues ||3^ 



B. Angular-momentum projection 



The practical techniques of the exact angular momentum projection were discussed many 
times in the literature [^,0. We review here the formalism in the form suitable for the 
description of the quadrupole degrees of freedom. 

Let us denote by \NI3'^) the particle- number-projected states ( |3.1| ) in the intrinsic frame 
defined by Eq. (|2.12| ), where /3 and 7 are the axial coordinates on the ^20^^22 plane 



C20 



[3 cos 7 



^2^22 = /? sin 7 . 



(3.9) 



According to standard prescriptions ||5,18.19,30.31], the angular-momentum-projected GCM 
state \nNIM) has the form 



\nNIM) = E [dPd-f ^"'""^f'^^ \NIMKp-f) 
J 1 -\- Oko 



where the angular-momentum-projected intrinsic state \NIMKj3'~f) is given by 

21 + 1 



(3.10) 



\nimkpj) = Jdn [dCk{q) + {-lYD^-Km 



RmN(3^) . 



(3.11) 



The rotation operator R{Q) depends on the three Euler angles, Q=(j), 6, ip. Due to the D2 
symmetry of the GCM states ( p.5|) the sum in Eq. ( p.lO|) is restricted to non-negative even 
values of K. Using states ( |3.10|) one obtains the GCM (Hill- Wheeler) equation in the form 
0J3|J3^ 



E 

K'>Q^ 



d(5'di 



n'KK'if^, 7; 7') - eMAP, 7; P', i) 



9nIKiP,7) = 



(3.12) 
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where Ti and are the GCM kernel matrices of the hamiltonian and of the norm, respec- 
tively, defined as 



2J + 1 



A 



KK' 



dn {Np-f\dRin)\Np'y) 



(3.13) 



X 



and Akk' stands for l/((l+5xo)(l+'^i^'o))- Since evaluation of kernels in Eq. (|3.13|) re- 
quires three dimensional integration, this is the most time consuming part of the numerical 
calculation. 

Because all indices K and K' are even, we may use the explicit form of the Wigner 
matrices, Dl^j^i{yL)=e^^^ S^j^i{9)e^^ ^ to rewrite ( p.l3|) in the following way 



21+1 



87r2 



-A 



KK' 



/■2-K 




/ 


b Fde sin 6 f 




Jo Jo 



2-K 



X 



dKK'{G) cos{K(j) + K'lP) + (-l)^rf^_x'(^) cos(K0 - K'lP) 



(3.14) 



By the same token the integration intervals over </> and can be reduced from (0, 27r) 
to (0,7r). Moreover, the Wigner functions d^j^j^,{9) depend then on a;=cos6' only, and the 

integral j d9 sin 9 can be rewritten as j dx. Furthermore, symmetry properties ||3^ of 

d^KK'i^) allow also to reduce integration over x to the interval (0,1). Finally, the integration 
domain in ( |3.14 ) is eight times smaller and the factor 2 in front is replaced by 16. 

Numerical integration over and ip Euler angles is done by means of the Tchebyshev 
method ||3^, i.e., values of the integrand in equally spaced points are summed up with equal 
weights. This method is exact when number of points is taken to be /max+l, where /max is 
the maximum spin of the good- angular- momentum components in the intrinsic state \N(3'^). 
Integration over x is done using the Gauss-Legendre method [^] applied to the (0,1) interval. 
Again, exact result is obtained by setting the number of points to Jmax+l- 



IV. GCM CALCULATIONS IN THE INTRINSIC FRAME 

The exact particle-number and angular-momentum projections allow us to calculate 
hamiltonian and norm kernels, 'H^KK'iP^I'^ P' ^l') J^KK'iPyl'^ P'^l')^ f*^^ arbitrary pairs 
of intrinsic deformations (/5,7) and (/5',7'), which define the GCM equation ( |3.12| ). By a 
discretization of the f3 and 7 variables this integral eigen-equation is transformed into a ma- 
trix eigen-equation, for which the kernels have to be calculated between all pairs of selected 
points {Pn,7n), n=l,. . .,Nd. 

In principle, the invariant measure sin 37|(i/3(i7 should be used in the GCM equation 
( p.l2| ) for quadrupole motion. However, since we discretize this integral equation, the inte- 
gration measure is immaterial. Moreover, we may then also include points for axial (7=0) 
and spherical (/3=0) shapes, which in the discretized GCM equation represent a certain 
volume of the phase space around points for which the invariant measure vanishes. 

For a given j, the main factor determining the total computing time is the number 
Nd of the mesh points (/3„,7„) on the /3-7 plane, because the number of projected kernels 
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to be calculated is equal hNd^N^ + 1). Each of these kernels requires calculating (/max+1' 



3 



overlaps for different Euler angles, Eq. ( 3.14 ). In Fig. ^we present the CPU times required to 



calculate one hundred of such overlaps as a function of the single-particle angular momentum 
j of the shell. Full circles denote results obtained on a 50MHz PC-486 scalar computer. This 
computing time grows roughly as a third-order polynomial in j (solid line), which is a much 
slower increase than what would be required for an exact solution. The times obtained 
on a vector supercomputer CRAY Y/MP EL98 are marked with full squares. They were 
multiplied by a factor of 20 before plotting them together with the PC results. The increase of 
time with j is here much slower, which is due to a better performance of the vector processor 
obtained for larger matrices. In fact, this time can be approximated by a polynomial of a 
second degree in j (solid line). 

In the present study we have performed the GCM calculations for j=31/2 and for j=15/2. 
For the latter value, the exact results can easily be obtained and were used to test the GCM 
method. 

As discussed in Sec. |II B| , the average value of angular momentum in the intrinsic states 
IN/Sj) increases with j3. Similarly, for large values of / the collective weight functions 
QniK^P-i'^) have large components at large [3. Choosing the mesh of points {(3n,ln) with 
smaller than a given value /3max is equivalent to restricting the variational space to the values 
of spins smaller than a certain maximum value /max- In calculations we used /3max=6 which 
allows to properly describe spins up to Jmax=18. 

Two meshes of points (/?n,7n) were used. The first one was composed of 10 points 
uniformly distributed in the sector /3</3max, 0<7<60° (Fig. |^, full circles) in such a way that 
overlaps of the neighboring points do not differ too much one from another. The average 
overlap turns out to be 0.73 for A^=10 particles in the j=31/2 shell and 0.80 for iV=8 and 
j=15/2. To test the influence of the number of points on the results, the second mesh 
comprising 19 points was constructed by adding one point in the middle of each triangle 
formed by three neighboring points of the 10-point lattice. 

In Fig. ^ the PES in the intrinsic frame of reference is presented for the pure quadrupole 
force. The upper part shows the results for iV=8 and j = 15/2, in which case the PES has 
a deformed 7-unstable valley. Two absolute minima are at /?~3.5 for 7=0° and 7=60° with 
a low saddle in between. The lower part of Fig. ^ shows the PES for A^=10 and j=31/2. 
In this case there appears a strongly pronounced oblate valley that extends towards large 
values of (3. 

The quadrupole moment of the intrinsic state can be calculated as an expectation value 
of the quadrupole operator in that state, 

Q,{P,^) = {NP^\Q^\NP^) . (4.1) 

The intrinsic frames of reference for the deformation tensor t2 and for the quadrupole tensor 
Q coincide, i.e., the only non-zero components of Q are 

Qo = Q cos 5 and V2Q2 = Qsm6 , (4.2) 



which defines the radial coordinates in the Q-6 plane, in analogy to Eq. (|3.9| ). When the 
mesh points of Fig. ^ are plotted in the Q-6 variables they appear to be distributed in a 
more uniform way than in the variables. On the other hand, the topology of the PES 
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is the same in both representations, with one important difference, namely, due to a finite 
dimension of the single-j-shell space, the quadrupole moment Q has its upper hmit Qmax- 
In Table ^ the values of Qmax are shown for several particle numbers in the j=15/2 and 
j=31/2 shells. 

For j =31/2 the intrinsic states \N(3'^) with the (3 coordinates beyond /5~10 all have 
quadrupole moments close to the maximum value Qmax- This means that there is no reason 
to extend the generating coordinates outside the /?max=10 circle, and in fact the value of 
/3max=6 used in the calculations is sufficient. 



A. GCM solutions for i=15/2 and iV=8 

In this section we present the GCM calculations for j=15/2 and N=8, where the exact 
solutions are available and can be used to test the properties of the GCM approach. The 
Hill- Wheeler equation ( |3.12| ) is solved by using the standard method of calculating the square 



root of the overlap matrix Mj^^' described in Ref. In what follows we compare the results 
obtained for 10 and 19 mesh points of intrinsic deformations, as defined in Sec. |rV|. 

As is well known, the GCM overlap matrix has the spectrum with the zero value as an 
accumulating point, cf. the discussion in Ref. [PSI. In principle this matrix is strictly positive 



definite but in practical calculations the negative eigenvalues may appear at the level of the 
numerical precision. This may be related to a finite precision of the diagonalization method, 
and especially to the numerical noise which may appear in calculating the elements of the 
overlap matrix. In Figs. |^(a) and |](b) we present the absolute values of the overlap matrix 
for J=0 and /=6, respectively, obtained with 19 mesh points. For /=0 the exact spectrum 
contains 7 states, see Fig. |^ or Table 0, and the overlap matrix has the same number of large 
eigenvalues separated by 10 orders of magnitude from the remaining 12 eigenvalues which 
represent the numerical noise. This illustrates a rare situation when the GCM generating 
functions exhaust the corresponding Hilbert space completely. 

A more typical spectrum of the overlap matrix appears for /=6, where the total number 
of eigenvalues equals to 76 (its dimension equals the product of the number of points (19) 
and of the number of K components (4), ii'=0,2,4, and 6). Even if this dimension is larger 
then the number of the exact /=6 states (31), the GCM generating functions do not exhaust 
the Hilbert space and no gap in the overlap spectrum occurs. In calculating the square root 
of the overlap one has to decide how many norm eigenvalues should be retained. Since the 
negative eigenvalues reflect the appearance of the numerical noise, their magnitude can be 
used as an indication of how much the results are perturbed by numerical uncertainties. 
Therefore, in the following we shall use and discuss a useful parameter denoted by /ccrit, 
which is equal to the number of norm eigenvalues larger than the absolute value of the 
smallest eigenvalue. In Fig. ^ we present a graphic illustration of how this parameter is 
defined. 

Let us denote by k the number of norm eigenvalues retained in solving the GCM equation. 
A choice of a proper value of k should be based on studying the stability of solutions when 
this parameter increases. In Fig. ^ we show the GCM energies for J=0 (full circles) obtained 
by using 10 (a) and 19 (b) mesh points. The exact energies are also shown (open squares) 
and placed between the results obtained with k=kcrit and fc=fccrit+l- It can be seen that 
the GCM calculation with k=7 give the exact spectrum with a very good accuracy. This is 
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due to the fact that that the complete Hilbert space is obtained in this case. For values of 
k=l,2, or 3 one obtains a fair description of 1,2, or 3 lowest states. On the other hand, the 
fourth state which first appears for k=4 is correctly split into two exact solutions only after 
all the 7 basis states are included. We may see this kind of effect in many other examples 
and one may interpret it as a manifestation of the coupling of collective and non- collective 
degrees of freedom. Within this interpretation, for k=A the GCM correctly singles out the 
fourth collective state, which in reality is fragmented among the forth and the fifth exact 
states. 

As soon as we include in the basis more states than the critical value k^rit, there appear 
in the spectrum spurious states which have nothing in common with the exact solutions. 
Such states appear in an unpredictable way at various energies and in realistic calculations 
cannot be distinguished from the real ones. In fact, they also pollute the wave functions 
of states which have energies seemingly stable when k increases beyond the value of fccrit- 
Therefore, one should always try to avoid using too large values of k and an identification 
of a highest acceptable value is essential [p^] . 

In Fig. ^ we present similar results for the 1=2 states. For both mesh sizes we obtain a 
good description of the four lowest states by using k=A. At higher energy one may identify 
two other structures which appear at k=5 and k=9, and for increasing k become fragmented. 
For the 10-point calculation, except of the four lowest states, none of the higher excited states 
is reproduced with /c</ccrit=12. On the other hand, for the 19-point mesh the exact spectrum 
is obtained for fc<A;crit=16, because the full Hilbert space is then again exhausted by the 
GCM generating functions. It is interesting to note that for 10 mesh points the fifth exact 
state appears only beyond the value of k=kcrit- This may illustrate an approximate character 
of our prescription to define the critical value fccrit- For 19 points the fifth state is stable 
already at k=k crit— 2. 

Based on this results we will suppose in what follows that using a relatively small number 
of mesh points of the intrinsic deformations (as compared to the number of exact states), 
one obtains a fair description of the collective subspace, and that only states with strong 
collective components are then present in the GCM spectra. By increasing the number of 
mesh points one includes noncollective configurations and the GCM yields collective states 
mixed with the noncollective ones. A separation of states into the collective and noncollective 
subspace is of course not a well defined procedure (both in the experiment and in the exactly 
solvable model discussed here). Which states are called collective is therefore a subject of 
a model interpretation. With increasing excitation energy the collectivity is gradually lost 
and the mixing with noncollective states increases. We may see, however, that the GCM 
can be used as a practical tool to select and decouple the collective space. 

In Fig. we compare the GCM spectra for spins up to /=6 with the exact results. 
Again, both the 10-point (a) and 19-point (b) results are shown. In the GCM calculations 
the cut-off parameters k were chosen to be equal to the critical values fccrit? which are listed 
in Table together with the corresponding numbers of exact states for a given spin. In 
all cases we have checked that the GCM states have correct average values of spin [P) and 
particle numbers, (A^) and (A^^), and that the average values of the hamiltonian {H) are 
equal to the GCM eigenenergies. These turn out to be useful criteria of identifying spurious 
states, because as soon as the spurious states start mixing with the physical ones the above 
average values become perturbed, even if the symmetries have been exactly restored when 
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calculating the GCM kernels. 

In the GCM spectra there is no 1^ states because they belong to representations of 
the point group D2 which are not included in the present calculations. Our choice of the 
generating functions ( p.5|) restricts them to the symmetric representation (+++) p4|, and 
therefore such is the symmetry of all GCM eigenstates, whereas the 1^ states belong to 
mixed-symmetry representations. In principle, some other high-spin many-particle states 
in the single-j shell may also belong to the mixed-symmetry representations and will not 
be accessible in our implementation of the GCM. This is not a case for the iV=8 and 1=2 
states in the j=15/2 shell discussed here because the GCM reproduces these states exactly. 
On the other hand, the GCM method gives much less odd-spin states then even-spin states 



as compared to the exact solutions. Fig. [T^. This can be understood by recalling that, 



for example, in the asymmetric rotor model |^ only (J— 1)/2 odd-spin states belong to the 
symmetric representation as compared to the J/2+1 even-spin states. One should also stress 
that the standard Bohr collective quadrupole model describes the symmetric states only. 

Comparing the results for 10 and 19 points. Fig. |ll], we may see that only three collective 
7=3 states are found in our GCM space in the former case. The lowest one clearly reproduces 
the lowest exact state while the first-excited state is probably fragmented among three 
excited exact solutions. For 19 points 7 states are found of which the lowest 4 correspond to 
the exact ones. For 1=5 only two collective states are found; the third excited state appears 
for 19 points at the correct energy, but it is absent in the 10 point results. Similar analysis 
shows that the four lowest 7=4 and 7=6 states can be considered as collective ones. 



The lowest GCM states with spins up to 7=18 are shown in Fig. |T^. Full circles denote 
states with positive signature (— 1)"^ and open circles with the negative signature. The states 
are connected by lines according to their stretched E2 matrix elements, see Sec. |1I A . In order 
to illustrate the convergence properties of the GCM we draw such lines only between states 
which have the GCM absolute energies precise to better than 3.5% as compared with the 
exact values. Again we see that the 10-point GCM calculation selects two lowest collective 
bands while in the 19-point results many higher excited, and not necessarily collective states 
are found. 



B. GCM solutions for j=31/2 



In Fig. [T^ we show the GCM energies for 7=0 states of A^=8 particles obtained for dif- 
ferent values of the cut-off parameter k. The vertical dashed lines are placed between the 
results obtained with k=kcnt and A;=fccrit+1 to show the values of the critical cut-off param- 



eters /ccrit determined as described in Sec. |1V A|. For the 10-point mesh in the intrinsic frame 



of reference. Fig. |T3|(a), we obtain the two lowest 7=0 states for k=l and 2, respectively. 
Only the first of them has a fairly correct energy already for k=l\ the energy of the second 
one decreases rapidly and stabilizes around k=Q. According to our method of identifying 
collective states, the third and probably the fourth collective 7=0 states appear at /c=3 and 
5, respectively, and then become strongly fragmented. This picture is confirmed by the 
19-point results presented in Fig. p!3|(b). 

A similar analysis of the 7=2 states. Fig. [13, shows that the lowest three of them have a 



collective nature while other three collective structures (appearing at A;=4, 8, and 11 in Fig. 
|Tl|(b)) become completely mixed and disappear in the dense non-collective spectrum. 
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Figs, ^(a-f) show the GCM spectra for particle numbers A^=4 to A^=14, respectively, 
obtained by using the 19-point mesh of intrinsic deformations. All presented GCM eigenen- 
ergies correspond to states fulfilling the criterion of stability with respect to the cut-off 
parameters fc<fccrit as well as have correct average values of (J^), (iV), (N^), and (H), as 
discussed in Sec. IIV A. 



The spectrum for N=4 can be compared with the exact results shown in Fig. ^ It is 
seen that the lowest GCM collective bands agree very well with the exact calculations. This 
is due to the fact that the dimension of the Hilbert space is here not too large and the GCM 
generating functions seem to exhaust it, at least for lower spins. The states for spins higher 
than 1=12 are reproduced less accurately and the bands gradually disappear, even if in the 
exact results they can be followed up to still higher spins. The bending up at higher spins 
of the n^=2 GCM bands illustrates the effect of decreasing possibility of the GCM states to 
exhaust the exact Hilbert space. 

For all particle numbers one may clearly distinguish at least five lowest collective bands. 
All of them are similar to those for N=4: and may accordingly be interpreted as the 7 
vibrations coupled to an axial rotor, see Sec. [11 A|. In all these cases one cannot see any 



candidate for low-energy (3 vibrations |23 



With an increasing particle number the GCM collective quadrupole bands loose their 
regularity. This effect can probably be attributed to a decreasing ability of the generating 
functions used in this study to exhaust the space of quadrupole collective excitations. But 
even for A^=14 the low-energy collective bands are visible. On the other hand, for the 7- 
unstable case corresponding to the half-filled shell (iV=16), we could not obtain any stable 
GCM solutions even by using the 19-point mesh of generating coordinates. 



V. SUMMARY AND CONCLUSIONS 

In the present paper we have discussed properties of states of particles moving in a 
single-j shell and interacting with the quadrupole-quadrupole force. This system exhibits 
collective as well as non-collective features and we have analyzed a possibility to select the 
collective subspace by using the generator coordinate method (GCM). 

The GCM generating functions have been constructed as coherent excitations built from 
the single-particle quadrupole operator. This type of building blocks has been proposed 
in the single-particle coherent excitation model (SCEM) which constitutes an interesting 
alternative to other models based on a building-block concept, e.g., the bosonic IBM model 



25| and the algebraic Ginocchio models ||4^. A collective model constructed in this way is 
based on purely fermionic objects and requires no closed algebraic properties. 

Depending on the number of different intrinsic-state configurations used in the GCM we 
obtain different precision of results as compared with exact solutions. When this number is 
relatively small, the GCM is shown to correctly select the quadrupole collective structures, 
while when this number is large all exact states are obtained with high precision. 

We have shown that in the single-j shell the low-lying quadrupole collective states can 
be interpreted in the frame of a standard geometrical collective model, either in its limit of 
an axial rotor coupled to 7-vibrations or in the limit of a 7-unstable motion. In the single-j 
shell we did not find collective bands of the beta-vibrational character, even for j as large 
as 31/2. 
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APPENDIX A: CANONICAL BASIS FOR TRANSITION MATRIX ELEMENTS 



Canonical basis for antisymmetric matrices has been introduced in Refs. pi| , ^ and, 
together with the Wick PD|fl and Onishi theorems, can be used to facilitate cal- 

culation of average values of fermion operators in quasiparticle states. Here we extend the 
concept of the canonical basis to non-diagonal (transition) matrix elements. 

We begin by considering the eigen-equation for the product of two complex antisymmetric 
matrices and C 



C^C) W = WD , (Al) 

where the columns of the matrix W are linearly independent vectors of the Jordan basis, 
and D is the Jordan block matrix with the eigenvalues Di at the main diagonal, the zeros 
or ones just above the main diagonal, and zeros elsewhere. In general, the eigenvalues Di 
are complex. 

Multiplying from the left-hand and right-hand sides the eigen-equation ( |Al| ) by and 
then transposing, we obtain that 

C (C+W-^^) = W-'^D^ , (A2) 

which multiplied by from the left-hand side gives 

(C+C) [C^W-'"') = [C+W-'"^) D'^ . (A3) 

We may now discuss several specific cases depending on the degeneracy of eigenvalues. 



1. Non-degenerate eigenvalues 

Suppose that the z-th eigenvalue is non-degenerate. Then the i-th column of the matrix 
C^W~^'^ is proportional to the i-th column of W , 

(C+iy-i^) =aWrm, (A4) 

where a is a proportionality constant. We keep the notation of indices m and m' denoting 
the single-particle states, although the results derived here do not depend on the assumption 
that these states belong to the single-j shell. Multiplying both sides of Eq. (^) by 
and summing over m we obtain 

a = (W-^C+W~^^)„ , (A5) 

which gives a;=0, because the right-hand-side is a diagonal matrix element of an antisym- 
metric matrix. From Eq. (K4) we now see that the i-th column of C^W~^ is equal to 



zero and hence from Eq. ( |A2| ) the corresponding eigenvalue is equal to zero. Therefore, 
there can be at most one non-degenerate eigenvalue and it must be equal to zero. All other 
eigenvalues must be at least pairwise degenerate. 
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2. Doubly degenerate eigenvalues with one eigenvector 



Suppose now that the matrix C^C has a doubly degenerate eigenvalue Di for which 
only one eigenvector exists, and let this eigenvector be equal to the i-th column of W . Let 
the ?-th column of W contains its partner in the 2x2 Jordan block corresponding to this 
doubly degenerate eigenvalue. From (|A1|) we then have 



and from 



C+C'W 
C+C'W 



c+w 



DiWml + W„ 



= Di (^C+W-^^ 



(A6a) 
(A6b) 

(A7a) 
(A7b) 



Comparing Eqs. ( |A6a|) and ( |A7a|) we see that the ?-th column of C^W ^'^ must be propor- 
tional to the i-th column of W , 



(A8a) 



while multiplying Eqs. (|A6b|) by a and subtracting from (|A7b|) we obtain that 



(A8b) 

where a and P are constants. However, summing Eq. ( |A8bD with W^l we obtain (3=0. 
Similarly summing Eqs. ( |A8a|) and ( |A8b|) with W~l and W^^, respectively, we obtain that 
a is simultaneously equal to both transposed matrix elements of an antisymmetric matrix 
, and hence a=0 too. Finally, both vectors at the left-hand sides of Eqs. ( |A8a| ) and ( |A8b| ) 
vanish, which is in contradiction with Eq. (|A^), where the D matrix has an non-diagonal 
matrix element equal to one. As a conclusion we rule out the possibility of one eigenvector 
existing for a doubly degenerate eigenvalue. Similar derivation can be repeated for higher 
degeneracies and we conclude that the matrix D is always diagonal, i.e., the number of 
eigenvectors of the non-hermitian matrix C^C equals to the number of dimensions. 



3. Doubly degenerate eigenvalue with two eigenvectors 

Now we have to again consider the case of the double degeneracy and analyze the struc- 
ture of the two corresponding eigenvectors. Let these eigenvectors be equal to the i-th and 
t columns of the matrix W. Similarly as in Sec. [A 1| we show that the i-th and t-th columns 



of the matrix C^W are also eigenvectors with the same eigenvalue and therefore must 
be linear combinations of the type: 



C+W-'^ . = aW^i + (3Wmi , (A9a) 

/ mi 

(C^W-^^) _ = + 5Wrni , (A9b) 
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where a, (3, 7, and 5 are constants. Summing Eqs. ( |A9a| ) and ( |A9b| ) with W~^'^^ and 
^~^'rrit ^6 obtain that a=0 and 5=0, respectively, while performing the summation with 
these factors reversed we obtain that /3=— 7. 

Similar procedure can be repeated for higher degeneracies with the same result that 
the eigenvectors of the matrix C^C are grouped in degenerate pairs, and that these pairs 
transform the matrix to its canonical form where it has the form of 2x2 antisymmetric 
blocks with zeros elsewhere. Therefore, this canonical form can be written as 



'1^"'"^) .. = SjCpj, , (AlO) 



where we use the convention of the state i being a partner of the state i, while i=i. The 
numbers Ci=ci are canonical matrix elements of the matrix while Sj=— are the phase 
factors, |sip=l, depending on the phase convention chosen for the canonical matrix elements 
Cj and on the phases of columns of the matrix W . 

Comparing Eqs. ( [Al[ ) and ( [A10| ) we see that the same matrix W simultaneously trans- 



forms the second antisymmetric matrix C to its canonical form 

(w^C'w).=s*d,5,,, (All) 

where the canonical matrix elements c- of C are related to the canonical matrix elements 
c* of C+ by 

cfc: = A . (A12) 

Since the numbers q depend on the normalization of columns of W , we may chose to work 
with the canonical basis for which Cj=l and c[=Di. However, a better choice is to use the 
normalization in which cf and c- are equal: 

c* = c: = /d;, (A13) 
with an arbitrary branch of the square root calculated for complex numbers Di. 



APPENDIX B: TRANSITION MATRIX ELEMENTS IN THE CANONICAL 

BASIS 

The Wick theorem allows to express non-diagonal matrix elements of fermion operators 
between the Thouless states ( |2.10| ) through three transition density matrices 



{C'\ \C) 

{C'\C) 
{C'\ \C) 
(C"|C) 

{C'\C) 



Pmm' — 



\l + C+C')-^C+C 

\i + c+c')-^c+ 

C\l + C+C')~^ 



(Bla) 
(Bib) 
(Blc) 



Using the canonical forms of the Thouless matrices and C", Eqs. (|A10|) and (|A11|) , we 
obtain the transition density matrices in the canonical basis: 
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«. - (M'-VIV)^. = ^ , (B2a) 



. (M'-.M'-^)^, = . (B2b) 



where the normahzation (|A13 ) was used. We see that in the canonical basis the pairing- 
tensor densities are related by R=—R'^, provided phases are chosen in such a way that the 
phase factors are real. 

We may now express transition matrix elements of fermion operators through their single- 
particle matrix elements in the canonical basis. For the one-body and the pair-transfer 
operators 

<5 = X! Qmrn'^m'^m' , (B3a) 
mm' 

0+ = E QLn'<^^' , (B3b) 
mm' 

= E ^mm'^m^m' ; (B3c) 
mm' 

-P = E Pmm'O'mO'm' , (B3d) 



mm' 



we have 

(C'\C) ~ ^mm' Pm'm — 1 -L £). ' 

\ I / mm' i ' i 

{C'\C) - ^, - E ' ^ ^ 

\ I / mm t ^ 

{C'\P+\C) _ V- p+ /+ _ V- PttsfVD't (r^Ac) 

(C'lC) l^^mm'l^m'm 2^ ^ , l^^CJ 

\ I / mm' t ' 

{C'\P\C) _ p _ PliSi\fDi (T\AA\ 

{C'\C) 2-/--"^-'- 1 + ^ ' ^^^^^ 

\ I / mm' I * 

where only diagonal or paired matrix elements of matrices 

Q' = W-^QW , (B5a) 

Q+ = W-^Q+W , (B5b) 

P+ = W-^P^W-^^ , (B5c) 

P' = W^PW , (B5d) 

have to be known. 

The use of the canonical basis also facilitates calculation of transition matrix elements 

of two-body operators. Moreover, the particle-number projection discussed in Sec. [Ill A| can 
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then be performed at a very low cost. This is so because the multiphcation of the Thouless 
matrix of the ket state by the gauge phase factor e^**^ does not change the canonical 
eigenvectors. Only the eigenvalues Di are multiplied by e^^"^. Here we only give expressions 
for particle-number-projected matrix elements of separable interactions: 



{C'\QQ+Pm\C) = E^ 



Q[MD,DjRlit,j) 



+ Q',MsjstJD,D,Rl{t,j 



{C'\P+PPn\C)= P;,P±s,s*JD,D^Rlit,j) 



+ 2P:^P;^D,D,Riit,j) 



(B6) 



(B7) 



where Rl{i,j) are the residues [|T5 



1 



2id 



1 + D,e2^<^ 



[B8) 



which can be calculated by using the discretized integrals discussed in Sec. |III A| . The overlap 
of Thouless states is given by 



(C"|C(0)) = det (l + e^''^C+C')^'^ = 11 + Ae''"^) 



(B9) 



i>0 



where notation z>0 means that only one state from every canonical pair is included in 
the product. 



APPENDIX C: MATRIX ELEMENTS BETWEEN TIME-EVEN STATES 

Suppose that the Thouless states are time-even, i.e., 

f\C) = \C) , (CI) 

where T is the antiunitary time-reversal operator, and suppose that the single-particle basis 
is closed with respect to the time-reversal, i.e., 

Ta^T^ = ^ U'^^iO^, , (C2a) 

m' 

famf^ = U^m'(^m' , (C2b) 
m' 

where [/ is a unitary antisymmetric matrix. Then the Thouless matrix defining the state 
\C) is time-even, namely, 

C+ = UC^U^ . (C3) 
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For every antisymmetric time-even matrix such as we may define the hermitian 
time-even matrix C, 

C = C+U* =C*U+ = UC^ = U^C , (C4) 

= C , C = UC*U+ , (C5) 

which contains the same information about the Thouless state \C) as the matrix C^. Results 
of Appendices ^ and ^ can now be significantly simplified. 

First, the matrix C~^C' defining all transitional matrix elements becomes quasihermitian, 

c+c = ca = ^1/2 (c^/^c'c'^'') C'^/^ , (C6) 

i.e., by a similarity transformation C*^/^ can be transformed into a hermitian matrix Z: 

Z = C^/^C'C'^'' . (C7) 

Therefore, all eigenvalues Di become real numbers. (Cases when C is singular can be 
considered separately) . 

Second, the eigenvectors of the hermitian matrix Z form a unitary matrix V, 

ZV = VD , (C8) 

and hence the eigenvectors of C~^C' (cf. Eq. ([Al|) ) are given by 

W = &I'^V . (C9) 



Third, the hermitian matrix in Eq. ( |C6|) is time-even and therefore has eigenvalues pair- 
wise degenerate due to the Kramers degeneracy. Therefore, only half of the eigenvectors 
have to be numerically calculated and the canonical pairs of Appendix ^ can be identified 
with the time-reversed pairs [l32 |. 
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TABLES 



TABLE L Absolute values of the reduced matrix elements of the quadrupole operator between 
the lowest 2'^ and 0^ states obtained in the exact calculation for j=15/2 and N=8. 



oj 0^ 0^^ o_l 0^ 

2f L921 0.017 0.003 

2^ L307 0.014 

2+ 1.231 0.034 

2t 0.086 0.790 0.208 

2^ 0.070 0.387 0.437 



TABLE IL Same as Table | for j=31/2 and N=4. 

ot ot ot o| 

2f 1.366 0.004 <0.001 <0.001 

2^ 0.099 0.109 0.001 <0.001 

2^ 0.005 1.091 0.018 0.002 

2| 0.001 0.181 0.206 0.015 

2+ <0.001 0.020 0.833 0.072 
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TABLE III. The decomposition of the A^=8 SCEM collective states (U) in the single-j shell 
for j = 15/2 into states of good angular momentum. The numbers of new states obtained for a given 
value Np as compared to those for Np—l are listed for every I. 
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TABLE IV. Maximum quadrupole moments 



__i Qma: 

15/2 8 1.6 

31/2 10 1.7 

31/2 12 2.0 

31/2 14 2.2 

31/2 16 2.4 
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TABLE V. Numbers of states iVexact of spin / in the single-j shell for j=15 /2 and A^=8 particles 

compared with the numbers /ccrit of eigenvalues of the overlap matrix retained for the 10 and 19 
point GCM calculations (see text). 

_[ iVexact fccrit (10 points) fccrit (19 points) 

7 7 7 

14- 

2 16 12 16 

3 13 3 7 

4 25 13 21 

5 21 6 12 

6 31 15 22 

7 26 8 14 

8 35 16 24 

9 29 10 16 

10 35 16 23 

11 29 11 15 

12 34 14 18 

13 27 10 11 

14 30 10 11 

15 23 6 6 

16 25 6 6 

17 19 3 3 

18 20 4 4 



26 



FIGURES 



FIG. 1. Exact yrast spectra for A^=8 and j=15/2 obtained with varying strength of pairing 
and quadrupole interactions. 

FIG. 2. Exact yrast spectra for for different particle numbers N and sheh angular momenta j. 

FIG. 3. Single-particle energies in the intrinsic frame for j=15/2 as functions of the quadruploe 
moment Q (the Nilsson diagram). Exchange term in the mean field is neglected. 

FIG. 4. Potential energy surfaces (PES) as functions of the /3-7 collective coordinates corre- 
sponding to j=15/2 and A^=8 (a), and to j=31/2 and A^=10 (b). Full dots denote 10 deformation 
points used in discretizing the GCM equation in the intrinsic frame. 

FIG. 5. Exact spectrum for A^=8 particles occupying the j=15/2 shell 

FIG. 6. Exact spectrum for N=4 particles occupying the j=31/2 shell 

FIG. 7. CPU times required to calculate one hundred of overlaps for different Euler angles as 
a function of the shell angular momentum j. Solid lines show polynomial fits to the PC-486 and 
CRAY results, 0.237j3 + 0.882j2 and 0.159j2 + 0.314j, respectively. 

FIG. 8. Absolute values of eigenvalues of the GCM overlap matrix for j=15/2 and A^=8. 
Results obtained with the 19-point mesh are shown for /=0 (a) and 1=6 (b). 

FIG. 9. GCM spectra of spin /=0 states obtained for different values of the cut-off parameter 
k. Results obtained with 10 (a) and 19 (b) mesh points of intrinsic deformations are shown for 
j=15/2 and N=8. 

FIG. 10. Same as Fig. | for the 1=2 states. 

FIG. 11. GCM spectra of low-spin spin /<6 states obtained with 10 (a) and 19 (b) mesh points 
of intrinsic deformations shown for j=15/2 and N=8. 

FIG. 12. GCM spectra of low-energy spin /<18 states obtained with 10 (a) and 19 (b) mesh 
points of intrinsic deformations shown for j=15/2 and A^=8. 

FIG. 13. Same as Fig. | for j=31/2 and iV=10. 
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FIG. 14. Same as Fig. for i=31/2 and N=W. 



FIG. 15. GCM spectra for j=31/2 obtained with 19 mesh points of intrinsic deformations for 
particle numbers between A^=4 and A^=14, (a) to (f), respectively. 
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